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INTRODUCTION 

A  common  approximation  for  the  speed  of  sound  has  been  to  use  the  change  in  pressure 
with  a  change  in  volume  along  the  Hugoniot  as  equal  to  the  change  in  pressure  with  volume  along 
the  isentrope.  This  approximation  has  been  used  successfully  but  is  not  strictly  correct.  In  order  to 
evaluate  the  accuracy,  the  following  calculations  comparing  the  partial  derivatives  of  pressure  to 
volume  along  the  isentrope  and  Hugoniot  were  conducted. 


COMPUTATIONS 


The  speed  of  sound  can  be  calculated  using  the  partial  derivative  of  pressure  with  specific 
volume  along  the  isentrope  by  rearranging  the  definition  of  the  speed  of  sound  in  its  most  common 
form:  the  partial  derivative  of  pressure  with  density.  From  the  definition  of  the  speed  of  sound: 


c  = 


(1) 


The  density  is  related  to  the  specific  volume  by: 


1 


P  = 


v 


(2) 


1 

dp  = - =dv 

v z 

Then: 


(3) 


c  = 


(4) 


Pressure  is  a  state  property  and  therefore  can  be  determined  from  a  function  of  two  other 
state  properties.  In  this  development,  specific  volume  and  energy  are  being  used. 


P  =  f  (v,  e ) 


(5) 


Then  by  the  chain  rule: 


Dividing  by  dv: 


dP  = 


dP 

dv 


dv  + 


dP 

de 


de 


dP 

dv 


dP 

dv 


dv 

e  dv 


dP  de 
de  v  dv 


(6) 


(7) 
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Which  is  the  general  solution  and  must  be  true  along  the  Hugoniot: 


dP 

dP 

dP 

de 

—  - 

H - 
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dv 
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dv 

e  de 

Vdv 

Or 


Developed  by  similar  methods 


dP 

dv 
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.) 


(8) 


(9) 

(10) 


(11) 


Equation  11  relates 


a »p 
dv 


to 


dP_ 

dv 


specifically  the  deviation  between  the  two  derivatives  given 


ov\s  ov\H 

by  1^1  (p-l  -^1  Y  Applying  definition  of  entropy  simplifies  the  computations. 

oe  hi  \ov\c  ov\u/ 


From  the  fundamental  thermodynamics  relation: 

de  =  Tds  —  Pdv 


(12) 


Along  an  isentrope,  ds=0 


(13) 


Further  simplifications  can  be  made  by  using  an  equation  for  the  energy  along  the  Hugoniot: 


1 

e  =  2?(vo  ~  v ) 


(14) 
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The  development  to  this  point  used  the  laws  of  thermodynamics  and  the  conservation 
equations.  As  such,  equation  12  will  hold  regardless  of  the  equations  of  state  used.  Utilization  of 
equations  of  state  can  further  simplify  this  equation. 

From  the  Mie-Gruneisen  equation  of  state: 


Y 

v 


=  a  constant 


(17) 


Then 


dP 

dv 


vP  dP 

yl.  ' 

= - —  +  -r— 

1 - -Wo  -  v) 

s  v  2  dv 

H 

v  2 

(18) 


Using  the  linear  equation  of  state,  us-up: 

us  =  co  +  up  0®) 

Which  when  combined  with  the  application  of  mass,  momentum,  and  energy,  the  pressure  of  a 
shockwave  is  calculated  as: 


Cq2Q 0  -  v) 

[v0  -  s(v 0  -  v)Y 


(20) 


dP  2  v0  +  s(l70  —  v) 

dv  H  C°  [v0  -  s(v0  -  v)]3 


(21) 


Using  equations  14  and  16,  a  comparison  can  be  made  between  the  variables.  Copper  was 
selected  as  a  basis  for  comparison,  and  the  following  variables  were  used  as  shown  in  table  1 
(ref.  1): 


Table  1 

Comparison  between  variables 


Density 

8930 

kg/mA3 

Ambient  speed  of  sound 

3940 

m/s 

Linear  us-up  coefficient  (s) 

1.489 

unitless 

Mie-Gruneisen  gamma, 
approximately  =  2s-1  (ref.  2) 

1.978 

For  comparison,  the  pressure  (fig.  1)  and  energy  (fig.  2)  were  plotted  against  volume  after  the  shock 
compaction. 
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Figure  1 

Copper  Hugoniot  pressure-specific  volume  plane 
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Figure  2 

Copper  Hugoniot  energy-specific  volume  plane 
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Figure  3  shows  a  comparison  between  the  rate  of  change  of  the  pressure  with  volume  along 
the  isentrope  and  the  Hugoniot.  These  derivatives  are  identical  at  the  initial  volume  and  diverge  a 
greater  amount  for  progressively  higher  compressions.  Figure  4  uses  the  data  from  figure  3  to 
compute  the  speed  of  sound  in  the  shocked  materials. 
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Comparison  between  rate  of  change  of  pressure  with  volume  along  isentrope  and  Hugoniot  in 
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Speed  of  sound  in  shocked  copper 
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It  has  been  proposed  to  use  the  compression  ratio  to  increase  the  speed  of  sound  for  the 
compressed  material.  Figure  5  compares  this  method  to  using  the  Hugoniot  for  computing  the  speed 
of  sound.  This  method  shows  a  greater  deviation  from  the  computed  speed  of  sound  when 
compared  to  the  utilization  of  the  Hugoniot. 


0.6  0.7  0.8  0.9  1 

VOLUME  RATIO 


Speed  of  Sound  from  Volume 
Ratio 

Speed  of  sound  from 
Hugoiot/lsentrope 


Figure  5 

Copper  ratio  of  speeds  of  sound 


CONCLUSIONS 

The  methods  investigated  for  approximating  the  speed  of  sound  (using  the  Hugoniot  or  a 
compression  ratio  factor)  have  a  high  degree  of  accuracy  at  low  pressures  and  correctly  converge  to 
the  ambient  speed  of  sound.  All  methods  correctly  predict  higher  speeds  of  sound  for  the  higher 
energy  shocked  states.  The  approximations  diverge  progressively  greater  at  higher  shock  states. 

A  number  of  equations  were  developed  that  relate  the  partial  derivatives  of  the  change  in 
pressure  to  the  volume  along  the  Hugoniot  to  the  isentrope.  These  equations  are  correct  only  for 
states  on  the  isentrope  (those  obtained  as  a  function  of  a  single  shock).  The  first  uses  only 
thermodynamics  and  the  conservation  equations  and,  as  such,  would  be  universally  applicable. 

The  second  assumes  an  equation  of  state  that  allows  for  the  computations  of  the  values  for  the 
speed  of  sound. 
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